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Abstract. We study relations among special values of zeta functions, invariants of toric 
varieties, and generalized Dedekind sums. In particular, we use invariants arising in the 
Todd class of a toric variety to give a new explicit formula for the values of the zeta function 
of a real quadratic field at nonpositive integers. We also express these invariants in terms of 
the generalized Dedekind sums studied previously by several authors. The paper includes 
conceptual proofs of these relations and explicit computations of the various zeta values and 
Dedekind sums involved. 
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1. Introduction 

In the present paper, we study relations among special values of zeta functions of real 
quadratic fields, properties of generalized Dedekind sums and Todd classes of toric varieties. 
The main theme of the paper is the use of toric geometry to explain in a conceptual way 
properties of the values of zeta functions and Dedekind sums, as well to provide explicit 
computations. 

Both toric varieties and zeta functions associate numerical invariants to cones in lattices; 
however, with different motivations and applications. Though we will focus on the case of 
two-dimensional cones in the present paper, we introduce notation and definitions that are 
valid for cones of arbitrary dimension. The reasons for this added generality is clarity, as 
well as preparation for the results of a subsequent publication. 

1.1. Zeta functions. We begin by reviewing the first source of numerical invariants of 
cones: the study of zeta functions. 

Given a number field K, its zeta function is defined (for Re{s) sufficiently large) by: 



^ 0(a)- 

where the summation is over all nonzero ideals, and Q is the norm. The above function admits 
a meromorphic continuation in C, with a simple pole at s = 1, and regular everywhere else. 
Lichtenbaum ||L| conjectured a specific behavior of the zeta function Ci^i at nonpositive 
integers related to the global arithmetic of the number field. In the special case of a totally 
real number field Lichtenbaum conjectured that the values of the zeta function at negative 
integers are rational numbers which involve the rank of the algebraic (or etale) iiT-theory 
of K. It turns out (see |^hlj| and ||Za2| , Section 2]) that for a totally real field K the zeta 



function can be decomposed as a sum C,{K^ s) = ({I, s) where the sum is over the (finite 
set of) narrow ideal classes. For each narrow ideal class /, there is a finite set of open 

cones in a rank [K : Q] lattice M/ in K such that ({I~^,s) = J2iCQ,(T° j{s) where Q is a 
multiple of the norm on K and where for an open cone r° we have 

Cq,As)= E Ql^s- 

The problem of calculating the zeta values (^q t-o(— n) for > for all triples (M, Q, r°) that 
arise from totally real fields has attracted a lot of attention by several authors. Klingen ||Ki | 



and Siegel [pil| , pi2|| (using analytic methods) proved that the values of the zeta functions of 
totally real fields at nonpositive integers are rational numbers and provided an algorithm for 
calculating them. Meanwhile, Shintani [^hl|] (using algebraic and combinatorial methods) 



gave an independent calculation of the zeta values Cq,t°(~'^) for an arbitrary Q which is a 
product of linear forms. Meyer ||Me|| and Zagier ||Za3|| gave another calculation of the zeta 



values Cq,t-° (—■"-) for rank two lattices. Related results have also been obtained by P. Cassou- 
Nogues, [P=M| p:N2| , [CTS^ . More recently, Hayes 0, Sczech |Sc| and Stevens 



have constructed PGLm(Q) cocycles which, among other things, provide a calculation of 
the zeta values (at nonpositive integers) of totally real number fields in terms of generalized 
Dedekind sums. 

We now specialize to the case of real quadratic fields. Given a narrow ideal class /, there 
is a cone o"/ (i.e., the convex hull of two rays from the origin) in a rank two lattice M O K, 
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such that: ({I ^, s) = (q,cti{s). Here Q is a multiple of the norm form and Cq,o- is given by 
the following definition: 

Definition 1.1. Let M be an two-dimensional lattice, and let V denote the real vector 
space M]R = M ® M. Let Q : — >■ M be a nonzero quadratic homogenous function i.e., a 
function satisfying Q{av) = a^Qiy) for a G M, f G and such that Q is not identically 0. 
We set 

wt{T, a) 



(1) cqAs)= J2 



QiaY 

where r is a rational two-dimensional cone in M and wt{T, ■) : M — >■ Q is the weight function 
defined by: 

{1 if a lies in the interior of r 
1/2 if a lies in the boundary of r, and a 7^ 
otherwise. 

We will assume throughout that Q is positive on all of r; that is, for all a G r, a 7^ 0, we 
have Q{a) > 0. The above zeta function, defined for Re{s) sufficiently large, ||Za4|| can be 



analytically continued to a meromorphic function on C, regular everywhere except at 1. 
In [ ^a4 |, Zagier showed that after possibly multiplying M by a totally positive number 



and Q by a nonzero rational number (which only multiplies values of the zeta function 
by a nonzero constant), every triple {M,Q,t) which arise from real quadratic field can be 
constructed by means of a finite sequence b = {bo, . . . , 6r-i) of integers greater than 1 and 
not all equal to 2, which will be denoted by (Mfe,Qb,rb). In addition, there are canonical 
vectors Aq, . . . ,Ar in (that depend on b) that subdivide the cone Tb = {Aq, Ar) (i.e., the 
cone whose extreme rays are Ao and A^) in Mb in r nonsingular cones {Ai,Ai+i). For a 
detailed discussion, see Section |3TT| . 

We will now give an explicit formula for the zeta values CQ,„Ti,{—'n) (for n > 0), which 
in turn enables one to calculate the values of the zeta function of a real quadratic field at 
nonpositive integers. While both the motivation and the proof involve concepts from the 
theory of toric geometry, the formula can be stated and understood without a knowledge of 
toric varieties. We will state the formula here, and in the next subsection, we will discuss 
concepts from toric geometry which are necessary for the proof and which lead to a conceptual 
understanding of the present formula. 

Let Am be defined by the power series: 

(3) = y: ^mh- 

thus we have: = {—l)"^B„Jm\ where is the m*^ Bernoulli number. (See also Defini- 
tion below.) Note that if m > 1 is odd, then A^ = 0. For d > 2 even, define homogeneous 
polynomials Pd{X, Y), Rd{X, Y) of degree d — 2 by: 

Pd{x,Y) = Yl A,A,x^-ir^-\ 

i+j=d,i,j>0 
vd—l I T^d—1 

R^{X,Y) = ^ ^ = X"'^ - X'^-^Y + ■■■ + Y''-\ 

A + y 
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where 62 = —I and 
We then have: 



1 for all even d > 2. 



Theorem 1. The values CQ^,T,,{—n) forn > are polynomials in hi with rational coefficients 
(symmetric under cyclic permutation of the hi) given explicitly as follows: 



(4) 



(5) 



g-Qb(xAi_i+yAi) 



^-Qb{xA^--i+yA^+i) 



The diamond symbol indicates the evaluation of a partial derivative at the origin. 
In particular, we obtain the formula of Zagier ||Zal| , Equation 3.3]: 



(6) 



1 

j=0 



1.2. Toric geometry. The second source of numerical invariants of cones comes from the 

or 



Fu 



theory of toric varieties. For a general reference on toric varieties, see fDa 

Founded in the 1970s, the subject of toric varieties provides a strong link between alge- 
braic geometry and the theory of convex bodies in a lattice. To each lattice polytope (the 
convex hull of a finite set of lattice points) is associated an algebraic variety with a natural 
torus action. This correspondence enables one to translate important properties and theo- 
rems about lattice polytopes into the language of algebraic geometry, and vice-versa. One 
important example of this is the very classical problem of counting the number of lattice 
points in a polytope. The early pioneers in the subject of toric varieties found that this 
problem, translated into algebraic geometry, becomes the problem of finding the Todd class 
of a toric variety. Much progress has been made over the past ten years on the Todd class 
problem, and this has led to a greater understanding of the lattice point counting question. 

One approach to computing the Todd class of a toric variety is the fundamental work 
of R. Morelli |[lVlo|| . He settled a question of Danilov by proving a local formula expressing 
the Todd class of a toric variety as a cycle. Another approach, introduced in by the first 
author in [0, |P2|| , is to express the Todd class as a polynomial in the torus-invariant cycles. 
Dedekind sums appear as coefficients in these polynomials, and this leads to lattice point 
formulas in terms of Dedekind sums, as well as new reciprocity laws for Dedekind sums [pi|| . 
Cappell and Shaneson ||CS|| subsequently announced an extension of the program of |pi|| in 
which they proposed formulas for the Todd class of a toric variety in all dimensions. In [p3 | 



it was shown that the polynomials of |pi| , |P2|| can be expressed nicely as the truncation of a 
certain power series whose coefficients were shown to be polynomial-time computable using 
an idea of Barvinok 



Ba 



A beautiful power series expression for the equivariant Todd class 
of a toric variety was given by Brion and Vergne in [pV2|| . Guillemin ||Gu|| also proved similar 
Todd class formulas from a symplectic geometry point of view. Furthermore, in |PV1|| , Brion 
and Vergne use the Todd power series of ||BV2|] to give a formula for summing any polynomial 
function over a polytope. 

The power series studied in [P3| and |pVl| , pV2|| are, in fact, identical (see Section 0) 



and play a central role in the present paper. A detailed discussion of the properties of these 
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power series, which we call the Todd power series of a cone, is contained in Section ^ In 
this section, we state two theorems about the Todd power series of a two-dimensional cone 
that we will need in our study of the zeta function. 

Given independent rays pi, . . . , p„ from the origin in an n-dimensional lattice N, the convex 
hull of the rays in the vector space = N^M. forms an n-dimensional cone a = (pi, . . . , p„). 
Cones of this type (that is, ones generated by linearly independent rays) are called simplicial. 
Let C"'{N) denote the set of n-dimensional simplicial cones of with ordered rays. There 
is then a canonical function 



t:C"(iV)^Q|xi,...,x„], 

invariant under lattice automorphisms, which associates to each cone a a power series t„ with 
rational coefficients, called the Todd power series of a. Several ways of characterizing this 
function are given in Section ^. These include an A^-additivity property (See Proposition 



271]) , an exponential sum over the cone (Proposition |2.3| ) and an explicit cyclotomic sum 
formula (Proposition \2.4\ .) 

In the case of a two-dimensional cone a, the coefficient of xy of to- (a;, y) was identified as a 



Dedekind sum [pi|]. Furthermore, in | |P1| | it was shown that reciprocity formulas for Dedekind 



sums follow from an iV-additivity formula for t. Zagier's higher-dimensional Dedekind sums 
ZaI] were later shown to appear as coefficients [B V2|. It is natural to conjecture that. 



in the two-dimensional case, all the coefficients of t are generalized Dedekind sums and 
that reciprocity properties of generalized Dedekind sums will be related to the A^-additivity 
formula of t. Indeed, this is the case: see Theorem 

We now present an explicit link between Todd power series and zeta functions. The 
following theorem expresses the values of zeta function at negative integers in terms of the 
Todd power series of a two-dimensional cone. Note that the idea of considering the Todd 
power series as a differential operator, and applying this to an integral over a shifted cone is 
not new: this was introduced in ||KP|| and developed further in [pVl|| . 



First we introduce some notation which is standard in the theory of toric varieties. If r is 
a cone in a lattice M, the dual cone f is a cone in the dual lattice = Hom(M, Z), defined 
by 

f = {v e N\ {v, u) >0 for all u G r}. 

The dual of an n-dimensional simplicial cone r = (pi, . . . , Pn) is generated by the rays Ui 
a normal to the n — 1-dimensional faces of r; hence f is also a simplicial cone. Given 
h = {hi, ... , hn) G M", we define the shifted cone T{h) to be the following cone in M^: 

'^W = ^ M]^\{ui, m) > —hi for alH = 1, . . . , n}. 

Here (and throughout) we have identified each ray ui with the primitive lattice point on that 
ray, that is, the nonzero lattice point on Ui closest to the origin. 

Given an n-dimensional simplicial cone r = (pi, . . . ,p„) in an n-dimensional lattice M, 
the multiplicity of r, denoted by mult(r), is defined to be the index in M of the sublattice 
Zpi + ■ ■ ■ + Zp„ (again identifying the rays pi with their primitive lattice points.) Thus the 
multiplicity of r is simply the volume of the parallelepiped formed by the vectors from the 
origin to the primitive lattice points on the rays of r. 
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Theorem 2. Let t be a two-dimensional cone of multiplicity q in a two-dimensional lattice 
M , and let a = f he the dual cone in N = Hom{M, Z). Then for all n > 0, we have: 

where all derivatives above are evaluated at hi = h2 = and 6nfi = 1 (resp. 0) if n = 
(resp. n 0). 

In this equation, (to-)d denotes the degree d part of the Todd power series thought of as 
an (infinite order) constant coefficients differential operator acting on the function h — > 
/^(/,) exp{-Q{u))du. 

The coefficients of for a two-dimensional cone may be expressed explicitly in terms of 
continued fractions. We now state this formula. 

Let (T be a two-dimensional cone in a lattice N. Then there exist a (unique) pair of 
relatively prime integers p, q with < p < q such that a is lattice equivalent (equivalent 

under the automorphism group of the lattice) to the cone cr(p =^ ((1, 0), (p, g)) in Z^. Such 
a cone will be called a cone of type (p, g). Let bi, hi, ki, Xi are defined in terms of the 
negative-regular continued fraction expansions: 

(7) - = [bi, . . . ,br-i], ^ = [bi,...,bi^i], Xi = -hiX + {qki-phi)y. 

P Ki 

(Throughout, we use such bracketed lists to denote finite negative-regular continued frac- 
tions.) We then have the following continued fraction expression for the degree d part {icr)d 
of the Todd power series to-. 

Theorem 3. For a of type {p, q) as above, and for d> 2 an even integer, we have: 

r r—1 

{ia)d{x, y) = qxy ^ Pd(-^i-i, Xi) + \dqxy ^ 6ii?d(Xi_i, X^+i) 

i=l i=l 

-XdixXt' + yXtl) + \5d,2qxy. 
If d >1 is odd, then {ta)d{x,y) = \Xd-iq'^~'^xy{x'^^'^ + y'^^'^). 

Remark 1.2. We have stated the formula above in the form we will need it for our study of 
zeta functions. However, from the toric geometry point of view, it is more natural to use the 



continued fraction expansion of instead. This formula will be given in Section 2.3 

q-p 

Remark 1.3. The formula for the Todd operator in Theorem |^ is reminiscent of formulas in 
quantum cohomology, even though we know of no conceptual explanation for this fact. 

1.3. Dedekind sums. Finally, we calculate the coefficients of the Todd power series of 
two-dimensional cones in terms of a particular generalization Sjj (defined below) of the 
classical Dedekind sums. For an excellent review of the properties of the classical Dedekind 
|R-G|| . Several generalizations of the classical Dedekind sums were studied by T. 



sum, see 



Apostol, Carlitz, C. Meyer, D. Solomon, and more recently and generally by U. Halbritter 
||AV| , [FIa| , |Me| , [5o| . These papers investigate the sums ( and several other generalizations of 
them, which we will not consider here) given in the following definition: 
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Definition 1.4. For relatively prime integers p, q (with q ^ 0), and for nonnegative integers 
i,j, we define the following generalized Dedekind sum: 

(8) = 

i.j. q q 

where for a real number x, we denote by (x) the (unique) number such that (x) G (a; + 
Z) n (0,1]. Here Ber^ denotes the rn}^ Bernoulli polynomial, defined by the power series 
Yl^=o^^'^rn{x)^ = Ti^; cind Ber^ denotes the restriction of the m*^ Bernoulli polynomial 

Ber^ to (0, 1], with the boundary condition Berm(l) =^ |(Berm(l) + Berm(0)) = i?m + 5m^i/2, 
where Bm is the m^^ Bernoulli number, defined by Bm = Berm(O). Q 

By definition, the sum Si^i{—p,q) coincides with the classical Dedekind sum s{p,q) (cf. 
[R-G| |.) An important property of generalized Dedekind sums is a reciprocity formula which 
leads to an evaluation in terms of negative-regular continued fractions. This reciprocity 
formula was most conveniently written by D. Solomon in terms of an additivity formula of 
a generating power series [^ , Theorem 3.3]. On the other hand, the Todd operator also 
satisfies an additivity property. We denote by fi,j{p,q) (for nonnegative integers i,j) the 
coefficient of x^y^ of the power series to-^^ (x, |/) {abbreviated hj t(j,g){x,y)). We then have: 

Theorem 4. Let p, g G Z be relatively prime with q 0. If i,j > 1, then we have: 

(9) kj{P,q) = q''''-\-iys.,{P,q). 

If i = 1 or j = 1, the above equation is true when the correction term 

is added to the right hand side. 

Corollary 1.5. Fixing r, for all nonnegative integers i,j and with the notation of (^, 
Sij{p, q) are polynomials in hi, 1/q. 



Remark 1.6. For i 



J 



1, the above theorem was obtained in FT , and was a motivation 



for the results of the present paper. 

1.4. Is toric geometry needed? Some natural questions arise, at this point: 

• Is toric geometry needed? 

• How do the statement and proof of Theorem |^ differ from the statement and proof of 
Zagier's formula ||Za4|| ? 

With respect to the first question. Theorem (an evaluation of zeta functions associated 
to real quadratic fields) is stated without reference to toric geometry, and a close examina- 
tion reveals that its proof is based on an analytic Lemma 2.11| and on the two dimensional 
analogue of the Euler-MacLaurin formula given by Proposition p.l3| . In addition, two di- 
mensional cones can be canonically subdivided into nonsingular cones, so in a sense the two 



^There seem to be two conventions for denoting Bernoulli numbers; one, which we follow here, is used often 
in number theory texts. Due to the facts that B2n+i = for n > 1, and that the sign of -8211 alternates, 
the other convention, which is often used in intersection theory, defines the n*'* Bernoulli number to be 

(-1)"+1S2„. 
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dimensional analogue of the Euler-MacLaurin formula can be obtained by the classical (one 
dimensional) one, as is used by Zagier [ Za4 |. Furthermore, number theory offers, for every 



m > 2, a canonical Eisenstein cocycle of PGLm(Q) ||Scl| , ^c2|| that expresses, among other 
things, the generalized Dedekind sums Sij in terms of negative-regular continued fraction 
expansions like the ones of Theorems ^ and ^, and the values (at nonpositive integers) of 
zeta functions of totally real fields (of degree m) in terms of generalized Dedekind sums. See 
also [Bol, Bt . 



On the other hand, toric geometry constructs for every simplicial cone a in an m-dimension- 
al lattice A^, a canonical Todd power series satisfying an A^-additivity property (see Propo- 
sition |2.1| below). The coefficients of the Todd power series are generalized Dedekind sums, 
and the power series itself is intimately related to the Euler-MacLaurin summation formula. 
As a result, for m = 2, we provide a toric geometry explanation of Theorems ^ and ^. 



With respect to the second question, Zagier ||Za4|| obtained similar formulas for the values 
of the zeta function of a quadratic number field at negative integers. Zagier used additivity 
in the lattice M, whereas we use additivity in the dual lattice A^. A^-additivity implies that 
these zeta values are a priori polynomials in bi and moreover these polynomials are 
directly defined in terms of the Todd operator and the quadratic form Q. Thus the present 
paper is an example of the rather subtle difference between A^-additivity and M-additivity. 
This difference can be seen by contrasting the additivity property of the zeta function (g^r 
with the additivity property satisfied by Todd operator. The modified Todd operator 5^ 
(introduced in Section is A^-additive under subdivisions of the cone a, i.e., it satisfies 
equation ( pH] ) below. On the other hand, the function Cq,-{^) '■ C^(M) — C is M-additive, 
i.e., it satisfies: 

(10) CQ,r(s)=CQ,n(^)+CQ,r,(s) 

whenever r is subdivided into cones ti and T2. This M-additivity of (q,.{s) is due to the 
particular choice of the weight function wt, and is purely set-theoretic. In contrast, the 
A^-additivity of is somewhat deeper, and has the advantage that only cones of maximal 
dimension are involved. 

1.5. Plan of the proof. In Section ^ we review well known properties of Todd power series 
and prove Theorems |^ and ^ In Section |^, we review the relation between the zeta functions 
that we consider and the zeta functions of real quadratic fields. We give a detailed construc- 
tion of zeta functions, and prove Theorem Finally, in Section ^ we review properties of 
generalized Dedekind sums and prove Theorem ^ 

1.6. Acknowledgments. We wish to thank M. Rosen for encouraging conversations during 
the academic year 1995-96. We also thank B. Sturmfels and M. Brion for their guidance. 
We especially wish to thank W. Fulton for enlightening and encouraging conversations since 
our early years of graduate studies. 

2. The Todd power series of a cone 

In this section we study properties of the Todd power series associated to a simplicial 
cone a. Section ^]l| provides an introduction and statements of several previously discovered 
formulas for the Todd power series. In Section |2.2| we make these formulas explicit for two- 



dimensional cones. Section |2.3| contains a proof of the explicit continued fraction formula for 



the Todd power series of a two-dimensional cone. Finally, in Section p.4| we prove Theorem 
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which hnks Todd power series with the problem of evaluating zeta functions an nonpositive 
integers. 



2.1. General properties of the Todd power series. Todd power series were studied in 
connection with the Todd class of a simplicial toric variety in |[P3|| . Independently, they were 
introduced in | |B V 2| | in the study of the equivariant Todd class of a simplicial toric variety. In 



addition, these power series appear in Brion and Vergne's formula for counting lattice points 
in a simple polytope ||B V 1[ | , which is an extension of Khovanskii and Pukhlikov's formula 



In these remarkable formulas, the power series in question are considered as differential 
operators which are applied to the volume of a deformed polytope. The result yields the 
number of lattice points in the polytope, or more generally, the sum of any polynomial 
function over the lattice points in the polytope. Below (Proposition 2.13|) , we give a version 



of this formula expressing the sum of certain functions over the lattice points contained in a 
simplicial cone. 

The Todd power series considered in the works cited above are also closely related to 
the fundamental work of R. Morelli on the Todd class of a toric variety. ||Mo|| . A precise 



connection is given in [p3| . Section 1.8]. 

We now state some of the properties of the Todd power series of a simplicial cone. Our 
purpose is twofold: we will need these properties in our application to zeta functions, and 
we wish to unite the approaches of the works cited above. Here, we follow the notation of 

El- 



Let a = (pi, . . . pn) be an ra-dimensional simplicial cone in an ra- dimensional lattice N . The 
Todd power series i„ of a is a power series with rational coefficients in variables Xi, . . . 
be corresponding to the rays of a. These power series, when evaluated at certain divisor 
classes, yield the Todd class of any simplicial toric variety [[F3|, Theorem 1]. 



To state the properties of to-, it will be useful to consider the following variant s of the 
power series t defined in [ P3[ ] by: 



1 



mult(cr)xi ■ ■ ■ X, 



-t(j {x\ , . . . , Xji) , 



which is a Laurent series in xi, . . . , 

The Todd power series and So- are characterized by the following proposition |P^, Theo- 
rem 2], which states that s is additive under subdivisions (after suitable coordinate changes), 
and gives the value of s on nonsingular cones. An n-dimensional cone is called nonsingular 
if it is generated by rays forming a basis of the lattice. It is well-known that any cone may 
be subdivided into nonsingular cones, and that such a subdivision determines a resolution 
of singularities of the corresponding toric variety (cf. |]Fu| , Section 2.6].) 

Proposition 2.1. If T is a simplicial subdivision of a then 

(11) So(X) = J2 5,(7"VX). 

7Gr(„) 

where the sum is taken over the n-dimensional cones of the subdivision T. Here X denotes 
the column vector (xi, . . . ,Xn)* and we have identified each cone (a and ^) with the n-by-n 
matrix whose columns are the coordinates of the rays of that cone. 
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For nonsingular cones, a, we have the following expression for s„: 

n ^ 



Remark 2.2. It follows immediately that for any cone a, 5o- is a rational function of the 

The Laurent series in xi, . . . , x„ may also be expressed as an exponential sum over the 
lattice points in the cone, in the spirit of the important earlier work of M. Brion |^r| . 



Proposition 2.3. Let a denote the dual cone in the lattice M = Hom{N, Z). We then have 



^-{{in,pi)xiA \-(m,pn)xn) 



The equality is one of rational functions in the 



Proof. As noted in |^|, since (m, pi) > for m G a fl M, the right hand side has a meaning 
in the completion of C[yi, . . . , ?/„] with respect to the ideal {yi, . . . , yn), where yi stands for 
g-a^i While it is not obvious, the right hand side is an rational function of the e^'. See 

By 



r|, p. 654]. The left hand side is a rational function of the e^' by Proposition 2A 



p. 655] and the second formula of Proposition these two rational functions are equal on 
nonsingular cones. As any cone can be subdivided into nonsingular cones, it suffices to verify 
that the right hand side satisfies the additivity formula of Proposition But this follows 
again from Brion's work. See the proposition of ||B^, p. 657], for example. Intuitively, this 
additivity can be seen from the fact that a sum of exponentials over a cone containing a 
straight line vanishes formally. □ 

The also have an explicit expression in terms of cyclotomic sums, due to Brion and 
Vergne. Following [[B V2|] , we introduce the following notation. Let ui,...,Un denote the 
primitive generators of the dual cone a. Thus we have {ui, pj) = if i ^ j, and {ui, pi) e 
Z, but does not necessarily equal 1. Let A^^. be the subgroup of A^ generated by the pi, 
i = 1, . . . , n, and let Go- = N/N„. Then G„ is an abelian group of order mult (a), which we 
denote by q. Define characters of G„ by 



Proposition 2.4. The Laurent series s^j coincides with Brion and Vergne 's formula express- 
ing their Todd differential operator. Namely, we have 



n 



Proof. By inspection, the right hand side is a rational function in the variables = e~^' 
which takes the value 1 when yi = for all i. Such rational functions embed into the comple- 
tion of the ring C[yi, ...,?/„] with respect to the ideal {yi, . . . , yn). To prove the proposition, 
it is enough to show that the right hand side and the right hand side of Proposition |2.3| define 
the same element of this completion. To do so, we expand the right hand side, getting: 

hi 



En 



ai{g)e~ 
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which becomes 

^ fci,...,A;„>0 g£Ga 

This last sum over Gg- is either q or 0, depending on whether a^^ ■ ■ ■ a,^" is the trivial character 
of Go- or not. Comparing the above with the right hand side of Proposition 2, we see that it 
suffices to show that a^^ ■ ■ ■ a^" is the trivial character of G^r if and only if there exists m E a 
such that {m, Pi) = ki for all i. This straightforward lattice calculation is omitted. □ 

Remark 2.5. The sum in the proposition above appears in the important work of Ricardo 
Diaz and Sinai Robins |PR| | . They use such sums to give an explicit formula for the number of 



lattice points in a simple polytope. Interestingly, their techniques, which come from Fourier 
analysis, are seemingly unrelated to the toric geometry discussed above. 

2.2. Properties of Todd power series of two-dimensional cones. In this section, we 
state properties of the power series for a two-dimensional cone a. 

It is not hard to see that if a is any two-dimensional cone, then there are relatively prime 
integers p, q such that a is lattice-equivalent to the cone 

a(p,,)''=^'((l,0),(p,g))cZ2. 

Here q is determined up to sign and p is determined modulo q. Thus we may arrange to 
have g > and < p < g. This gives a complete classification of two-dimensional cones up 
to lattice isomorphism. In discussing Todd power series we will abbreviate ia^j,^) by i(p,q). 
The explicit cyclotomic formula of Proposition p.4| may be written as: 



Proposition 2.6. The Todd power series of a two-dimensional cone is given by 

xy 



(1 - cj^Pe-^)(l - ooe-y) 



Proof. With the coordinates above and the notation of Proposition [^.4| , we have Ui = (g, —p), 
U2 = (0, 1), and consists of the lattice points (0, /c), /c = 0, . . . , g — 1. The desired equation 
now follows directly from Proposition \2.4\ . □ 



In the two-dimensional case, the additivity formula of Proposition |2.1| can be expressed 
as an explicit reciprocity law. This, and a periodicity relation for s are contained in the 
following theorem. 

Proposition 2.7. Let p and q be relatively prime positive integers. Then 

pi g 1 

(12) S(p,g)(a; - -y, -y) +S(q,p)(y - -X, -x) = S(o,i)(a;, 

g g p p 

(13) S(p + g,g)(X,?/) = S(p,g)(X,?/) 

Proof. The quadrant ((1, 0), (0, 1)) may be subdivided into cones 71 = ((1, 0), (p, g)) and 
72 = ((0, 1), (p, g)). The cone 71 is of type (p, g), and 72 is of type (g,p). Applying the 
additivity formula (Proposition |2.1|) to this subdivision yields the first equation above. 

The second equation follows from the fact that the cones ((1, 0), (p, g)) and ((1, 0), (p+g, g)) 
are lattice isomorphic. □ 

Let s'^^ (resp. 3°'^'^) denote the part of s of even (resp. odd) total degree. 
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Corollary 2.8. Given a two-dimensional lattice N, the function s : C'^{N) Q{{x,y)) 
(where Q((x, y)) is the function field of the power series ring Q[x, y}) is uniquely determined 
by properties (|1^) and ([I3| ) and its initial condition S(o,i)(a;, ?/) = 1/((1 — e~^)(l — e~^)). 
Since equations ([T^ ) and (|13|) are homogenous with respect to the degrees of x and y, it 
follows that 5^^ (resp. s°'^'^j satisfies (JT^) and (JT^) with initial condition s^^ -^-^ (resp. s"^q\-^)- 

Remark 2.9. Corollary p.8| has a converse. iFiom. equations (|12]) and ([T3|) it follows that 
S(o,i) satisfies the following relations: 

5(0,1) (a;, y) = 3(0,1) (a; -y,y)+ S(o,i)(2/ - x, x) and 5(o,i)(x, y) = S(o,i)(i/, x) 

Conversely, one can show that given any element g{x,y) G Q((x, y)) satisfying: 

gix, y) = gix -y,y) + g{y - x, x) and g{x, y) = g{y, x) 

there is a unique function g : C'^{N) — > Q((x, y)) so that S(o,i) = 9- 

2.3. Continued fraction expansion for the Todd power series. In this section, we 
prove Theorem 1^, which expresses the coefficients in the Todd power series of a two-dimensional 
cone in terms of continued fractions. Before doing so, we first formulate an equivalent ver- 
sion which is more natural from the point of view of toric varieties. The continued fraction 
expansion in this second version of the formula corresponds directly to a desingularization 
of the cone. 

Given relatively prime integers p, q with p > and < p < g, let a,, 7,, 5,, Lj be defined 
in terms of the negative continued fraction expansions: 

— ^ = [ai, . . . ,a^_i], = [ai,...,ai_i], Li='jiX + {q6i + {p-q)'yi)y. 

q-p Oi 

We then have the following continued fraction expression for the degree d part (tro)d of 
the Todd power series of a two-dimensional cone w. 

Theorem 5. For w a cone of type {p,q) as above, and for d > 2 an even integer, we have: 



itua)dix, y) = -qxy ^ Pd{Li-u U) - Xdqxy ^ aiRd{Li_i, L^+i) + Xd{xLf ^ - yL'^^l) 

i=l i=l 

Ifd> 1 is odd, then {i^)d{x,y) = \\d-iq'^~'^xy{x'^-'^ + y'^''^). 

Proof. It will be convenient to choose coordinates so that w = ((0, —1), {q,q — p)) (which 
is easily seen to be lattice equivalent to the cone ((1,0), {p,q))). We now subdivide zu into 
nonsingular cones. It is well-known that for two-dimensional cones this can be done in 
a canonical way, and that the resulting subdivision has an explicit expression in terms of 
continued fractions ||Fu| , Section 2.6]. In our coordinate system, the rays of this nonsingular 
subdivision of m are given by 

/?o = (0,-1) 
Pi = (1,0) 
p2 = (ai,l) 



Ps = 



{q,q-p). 
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Thus we have + = aj/?j, which imphes that Pi = (7^, 6i). 

The cone w is subdivided into cones Wi = , i = 1, . . . s. The A^-additivity formula 

of Proposition |2?^ imphes that: 

s 

where we again we have identified the 2-dimensional cone with the 2-by-2 matrix whose 
columns are the primitive generators of w. One easily sees that this becomes 

s 

s^ix,y) = ^s^^(-Li_i,Li) 

i=l 

Rewriting the equation in terms of the yields 

y) = -qxy > 

1=1 

Since every Wi is nonsingular, we have 

i^SX,Y) = giX)giY) 

where g{z) = z/ (1 — e~^). Now consider the degree d part of the above. We will assume d > 2 
is even and leave the other (easier) cases to the reader. The degree d part of g{Li_i)g{Li) 
may be written as: 

—Li-iLiPii{—Li_i, Li) + \ii{Lf_^ + Lf). 

Summing the first term above yields the first term in the equation of the theorem. So it 
suffices to examine the remaining term: 

I— 1 

Using the relation Lj_i + Lj+i = ajLj, the sum above may be rewritten as 



a., 



i=l 



Li-i + -Lj+i \ Lq L 



Keeping in mind that Lq = ~qy and Lg = qx, Theorem ^follows easily. □ 
Corollary 2.10. For a two-dimensional cone a = {pi-, P2) of multiplicity q we have: 

(14) t.(/ii, h^) - \iip,){qh,)h2 - \ii,,){qh2)h, = Cih, h2) - ^ 

where is the even total degree part of the power series t^. 

Proof. It follows immediately from (and in fact is equivalent to) the formula for the odd part 
of the Todd power series tg- given by Theorem |^. □ 

We now prove Theorem ^. Let a be a two-dimensional cone of type {p, q) in a lattice A^, 
and let ai,hi, ki and Xj be as in Theorem ^. The dual cone a in the dual lattice M is easily 
seen to have type {—p, q) and so we may choose coordinates in M so that a = ((0, —1), {q,p)) 
in Z^. Furthermore, the negative- regular continued fraction expansion of q/p corresponds 
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naturally to the desingularization of the dual cone a. Explicitly, the desingularization of & 
is given by the subdivision 

Po =(0,-1) 
Pi =(1,0) 

P2 =(&1,1) 

Pr = (g,p)- 



One has p^+i + pj_i = hipi and thus pi = {hi, ki). Applying Theorem ^ to d" expresses i{-p,q) 

le re. 

qxy 



in terms of the 6j. However t(-p,g) is related to i(p.q) via the relation 



1 - e-iy ■ 



In this way, we obtain an expression for t(p^g) in terms of the 6j, which concludes the proof 
of Theorem ^ 

2.4. Zeta function values in terms of Todd power series. In this section, we prove 
Theorem ^ which expresses values of the zeta function of a two-dimensional cone in terms 
of the Todd power series. Three ingredients are involved in the proof of this theorem: an 
asymptotic series formula (Proposition p.ll| ), a polytope summation formula (Proposition 
2.12| ) and a cone summation formula (Proposition |2.13|) . We begin with the first ingredient: 



Lemma 2.11. ||Za4| , Proposition 2] Let 0(s) = X1a>o '^^-^^^ ^ Dirichlet series where {A} 
is a sequence of positive real numbers converging to infinity. Let E{t) = Ylx>o^>^^^^^ 
corresponding exponential series. Assume that E{t) has the following asymptotic expansion 
as t ^ 0: 



(15) E{t) ~ Yl 



n=-l 

Then it follows that 

• (f){s) can be extended to a meromorphic function on C. 

• 0(s) has a simple pole at s = 1, and no other poles. 

• The values of (p o-t nonpositive integers are given by: (pi^—n) = (— l)"n!cji. 

We now present our second ingredient, a polytope summation formula. This proposition 
is a variant of the lattice point formula of [[BV2|| . A polytope of dimension n is called simple 
if each of its vertices lies on exactly n facets {n — 1-dimensional faces) of the polytope. 

Proposition 2.12. Let N be an n-dimensional lattice, P a simple lattice polytope in M and 
S its associated fan in N. For every analytic function : Mr M, we have the following 
asymptotic expansion as t — > 0; 

(16) V (j){ta) ~ ts ( ^ ) o / (P{tu)du 

where ts is the Todd power series of ||BV2| , Definition 10]. 
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Proof. First of all, the meaning of the right hand side is as follows: we consider the degree 
k Taylor expansion (p = (pk + Rk of 0, where (pk is a polynomial in Mr of degree k and Rk is 
the remainder satisfying lima^o |a|'^-RA;(fl) = 0. It follows that /p^^-, Pk{tu)du is a polynomial 
in t and h of degree k (with respect to t) and that /p^^^ Rk(tu)du = o{t^) at t = (with the 
notation that f(t) = o{t^) if and only if limi_o f(t)t~^ = 0). Thus, 

On the other hand, 

(f){ta) = Mta)+ Yl ^-^^^"^ 

a€PnM aePnM a€PnM 



aePnM 



Brion-Vergne |[BV2| , theorem 11] prove that for every poljTiomial function (such as (pk) on 
Mm we have: 



(pk{ta) = ts (St) o [ (pk{tu)du 



aePnM 

which concludes the proof. □ 

We call a function : — > M rapidly decreasing if it is analytic and for every constant 
coefficients differential operator D, and every subset I of [n] = {1, . . . ,n}, the restriction 
D{(p)\i obtained by setting a;^ = for i ^ / is in L^(M:!^). Examples of rapidly decreasing 
functions can be obtained by setting (p = Pexp(Q) where P is a polynomial on M" and 
Q : M" ^ R is totally positive, i.e., its restriction to M.^_^ takes positive values for every 
subset / of [n]. 

Proposition 2.13. Let N be an n- dimensional lattice. For every a G C"(A^), and every 
rapidly decreasing analytic function (p : Mk M, we have the following asymptotic expansion 
as t ^ 0: 

(17) V cPita) ~ (^) o [ cP{tu)du 

Proof. First of all, the right hand side of the above equation has the following meaning: 
consider the decomposition t„ = J2k ^fc '^^ power series t^, where ta-^k is a homogenous 
polynomial of degree k. A change of variables v = tu implies that 

ill) ^ L ^ (^) ' L ^ (^) ' L 

is a multiple of t'^"", thus the right hand side is defined to be the Laurent power series in t 
given by 



k=0 ^ ^ •^'^ 



(p{v)dv. 

(h) 
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For the proof of the proposition, truncate in some way the cone a in M to obtain a sim- 
ple convex polytope P, with associated fan S. Since Ur->o'"-P = cr, using the convergence 
properties of 0, we obtain as r — oo: 

0(ta) = hm (pita) 

~ limts (^-^j^ ^ [ (t){tu)du 

d_ 
dh 



J{rP){h) 

hmts o f <P{v)dv/e 

\OhJ J(^rtP){th) 

hmVts,fct'-"f^)o / 4>{v)dv 

^ V ^^^^ JirtP)ih) 



l{rtP)(h) 

o / (j){v)dv, 
Jp{h) 

which concludes the proof. □ 

In case of a two-dimensional lattice and a rapidly decreasing function : R, 
using the weight function wt of equation (j^) and inclusion-exclusion, we obtain the following 

Corollary 2.14. For a two-dimensional cone o = {pi, P2) of multiplicity q in N, we have 
the asymptotic expansion as t 0: 

(18) ^£„,...),,a)^{c(|^.^) -|^}o/^^_,W... 

where the right hand side lies in the formal power series ring t~^R[t^]. 

Proof of Theorem ^. Recall that r is a cone in M, a is its dual in and Q is homogenous 
quadratic, totally positive on r; thus e"'^ is rapidly decreasing on Mk. Corollary |2. 14] implies 
that the generating function 

(19) ZQAt)= ^tir,a)e-''^^''^ = wt{T,a)e-'^^'''"^\ 

satisfies the hypothesis of Lemma 2.11 , which in turn yields Theorem 0. □ 

Remark 2.15. Notice that the above proof of Theorem used crucially the fact that Bi = — ^ 
and the definition of the weight function wt. If we had weighted the sum defining the zeta 
function in any other way, the resulting variation of Theorem |^ would not hold. 

3. Zeta functions of number fields 

3.1. Parametrizing triples (M„, Q^, r^). In this section, we review in detail the construc- 
tion of the triple {Mj,, Qh, Tf,) mentioned in Section ^TT] that is associated to a sequence 
b = {bo, . . . , br~i) of integers greater than 1 and not all equal to 2. Our notation and con- 
struction is borrowed from [[Za4|| . Given a sequence b as above, we extend it to a sequence of 
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integers parametrized by the integers by defining bk = bkmodr- Furthermore, for an integer 
fc, we define 

Wk = [bk, . . . , 6fc+r-i] = h 



Jk+l 



bk+2 — ■ ■ ■ 

where [fefc, . . . , 6^+^-1] denotes the infinite continued fraction with period r. Note that Wk 
Wk+r for all integers k, and that, by definition, 

1 

Wo = bo 

61 



1 

■ ■ br-l 

Wo 



from which it follows that wo satisfies a quadratic equation Abw'^ + BbW + Cb = where 
Ab, Bb, Cb are (for a fixed r) polynomials in 6j with integer coefficients.^ Let Db = B^ — AAbCb 
be the discriminant. Since bi > 1, it follows that Db > 0, and that the roots of the quadratic 
equation are Wq = {—Bb + y/Db) / i'^Ab) and w'q = {—Bb — y/Db) / i'^Ab) . We thus have that 
Wq > 1 > w'q. We define the complex numbers Aq = 1, Ak-i = AkWk for k E Z. Since 
Wk = bk — it follows that 

(20) Ak^i + Ak+i = bkAk 

for all integers k. Let us define Mb = Zwo+Z to be the rank two lattice with the oriented basis 
{wq, 1}. Note that Mb is a rank two lattice in the real quadratic field Kb = Q(\/^). Due to 
the fact that A_i = wq and the recursion relation (PO]), it follows that Mb = ZAk + ZAk+i 
for all integers k, and thus eMb = Mb where e = wq . . . Wr-i- Thus Ub = {e^\m G Z} acts on 
Mb by multiplication, and the zeta function of (the inverse of) the narrow ideal class of Mb 
satisfies 



aMb,s)=N{Mb 



^ Qb{ay 



a>0 

where in the above summation we exclude 0, and N{Mb) is a constant, defined in |[Za4 
which is not so crucial for our purposes. We are therefore interested to describe the quotient 
Mb/Ub- The action of Ub on Mb implies that eAk = Ak-r- The matrix of multiplication by e 

on Mb, in terms of the basis {Pi, (32} == {^-1,^0} can be calculated as follows. Recall first 
that in terms of the above basis of Mb, we have that: 

Pk 

(21) Ak = -pkA_i + QkAo where — = [60, • • • , fefc-i] 

Qk 

which implies that 

(3i\ ^ M,_A ^ f-Pr-l -Pr\ fPl 

h) \A, ) \ ) \(32 

from which it follows that the matrix of e is given by 

|c.c.^ I Qr Pr \ _ f hq - p Q \ _ [^0 "A fK-l -l" 

' ' ' -Qr-i -Pr-1 ~ [-boP' + ipp' - l)/q -P' ~[l i • • • I 1 



^not to be confused with Aq, . . . , Ar introduced below, despite the rather poor choice of notation. 



18 



STAVROS GAROUFALIDIS AND JAMES E. POMMERSHEIM 



where 

(23) ^'='[6i,...,6,_i], 

(with the understanding that g = l,p = if r = 1), p' = numerator . . . , br-2], which is 
a multiphcative inverse of p mod q (with the understanding that p' = (resp. 1) if r = 1 
(resp. 2)). 

Lemma 3.1. The orbit space Mh/Ub has as a fundamental domain a semiopen two-dimensional 
simplicial cone whose closure is the cone {Ao,Ar) in Ki,. Furthermore, the (closed) cone 
Tb = {Aq, Ar) in Mb is of type {—p, q) and the dual cone ab in Nb is of type {p, q) where 
< p < q as in equation (^) above. 

Proof. /^From Figure |I], it follows that a fundamental domain of Mb/Ub is a semiopen two- 
dimensional cone whose closure is the cone Tb = {Ar,Ao) in Mb. The cone (^^,^0) is 
canonically subdivided into nonsingular cones {Ai^i,Ai) (for i = 0, ... ,r — 1), and using 
equation (|20|) , it follows that Tb is of type (ci, q) where < ci < g and = [61, ... , br-i]. 
Therefore the dual cone ab in is of type {p,q) where p = q — ci. Thus {p,q) satisfies 
equation (|23|) . □ 



A2 
^A, 

Ao = l = (32 
^A^i =wo = Pi 

Figure 1. A fundamental domain of Mb/Ub- 

Since Mb C Q{^/D^) is a rank two lattice of a real quadratic field, we can define a 
homogenous function Qb : Mr — M by 

Qb{xwo + y) = CbX^ - Bbxy + Aby"^ 

This agrees with the (normalized) norm function of the real quadratic field Qly/Db) consid- 
ered by Zagier | [^a4| , p 138]. We close this section with the following: 

Proposition 3.2. Let Z,m G Z with I ^ m. Then the coefficients of x'^,xy and y^ of the 
quadratic form Qb{xAi + yAm) are given by polynomials in bi with integer coefficients. 

Proof. We first claim that in terms of the basis {A^, Aq} we have: 

(24) Qb{xAr + yAo) = qQe{x, y) 

where Qe{x, y) = x"^ + Qxy + y"^ and 6 = bQq — p — p' . 

The above claim follows from the fact that Ab, Bb, Cb can be calculated (using their defi- 
nition) in terms of the bi as follows: 

Ab = q, Bb = -boq + p - p' , and Cb = qr-i = bop' + (1 - pp')/q, 

together with a change of variables formula from {y4_i, Aq} to {y4r,v4o} given by equation 
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We then claim that if l,m are distinct integers, then 
(25) QbixAi + yA„,) = -Qe{Hix + HmV, Kix + KmV) 

q 

where ^ = [61, . . . , HiX + Kiy = Xi = -hiX + (qki - phi)y. 

This follows from equation using the change of variables given in equation (pi]). 

We now claim that the coefficients of x'^,xy and y"^ in equation p5|) are polynomials in 
the bi with integer coefficients. 

Indeed, observe that hi, ki are polynomials in hi with integer coefficients. Since Hi = 
—hi,Ki = qki — phi = —phi mod q, it follows that: 

Qe{Hix + H^y, Kix + K^^y) = {hix + h^yf + QQnx + h^y){phix + ph^y) 

+ {phix + phrayY q 
= {hix + hmVy {I + Qp + p'^) mod q. 

Since 9 = qbo — p — p', and pp' = 1 mod q, it follows that 1 + Qp + = mod q, which 
finishes the proof of the proposition. □ 



Remark 3.3. The notation of Theorem |^ (and its proof, in Section p.3| ) corresponds to the 
notation of this section. Indeed, the cone a^, in is of type {p, q) where p, q are given in 
terms of equation (|23|). The rays pi (for i = 1, . . .r — 1) that desingularize the cone in 
M are given by pi = Ai. In addition, the hi,ki,Xi that appear in Theorem ^ match the 
notation Proposition p^ . 

3.2. Proof of Theorem g. 

Proof. The main idea is to use Theorem |^ which calculates the zeta values in terms of the 
Todd operator of cJb, and Theorem ^ which expresses the Todd operator in terms of the bi. 
The expression that we obtain for the zeta values differ from the one of equation (^ by an 
error term, which vanishes identically, as one can show by an explicit calculation. 

Now, for the details, we follow the notation of Section We begin by calculating the 
integral f , , e~^^'^^du. 

Using the parametrization — > Mr given by: (^1,^2) — » WiAr + W2A0, it follows that 
the preimage of Tb{x,y) in is given by {{wi,W2)\wi > —x/q,W2 > —y/q}- Thus we have: 
f 1*00 1*00 

JTt{x,y) Jwi=-x/q Jw2=-y/q 

Differentiating, we get 

g-QftW^y — l^-Qbi-xAr/q-yAo/q) _ lQ-l/qQe(x,y) 

{x,y) 

Using the following elementary identity: 
f d dV f d Oy _ _ f dV f dy 
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(and temporarily abbreviating ^ by x) we obtain that 

Jrhlx.v) 



'Tb{x,y) 

dx J \dy 

'd_y (d_ 

dx J \dy 



— \ f — \ ^ ^-l/qQe(Hix+Hmy,Kix+Kmy) 



Y J_ V o e-Q''(-^'+J^^'") , as well as 



{xX^^-^'+yX^^t')o I 

J tJx 



Q^^^ Jo \\dx dy J \dx dy 

The above, together with Theorem ^ implies that: 

CQb,n{-n) = {-irn\ <{ > ^2^+2 ( ^, ^ ) o e-«^(-^-^+^^') 



where 



i^;2„+2(&) = A2„+2{-&oi?2„+2(^,^)e-^/^(^^^+^^^+^^^) 



1 rlfd d^^' fd ,5^'"+' 



and 

L = qp'bo + l-pp', M = qboBb + 2{pp' - 1), = qpbo + 1 - pp'. 



ix 



Fixing r, the length of the sequence (60, • • • , K-i), Lemma below can be used to express 
E2n+2{b) as a polynomial in l/q,p,p',bo. An explicit but lengthy calculation implies that 
E2n+2{b) = for any b. We have thus succeeded in expressing CQb,(Tb{.~n) as a polynomial in 
bi with rational coefficients, symmetric under cyclic permutation, (since Qb{xAi + yAm) is 
symmetric under cyclic permutation). This concludes the proof of Theorem 0. Equation (^) 
follows easily using Ai = 1/2, A2 = 1/12. □ 



Lemma 3.4. [[Za4|| For every i,j>0 with i + j even, we have: 
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Furthermore, for every n> 1 we have: 



2c" r!(2n-l-2r)! 



— n \ t^X / 2c^ 



4. DEDEKIND SUMS IN TERMS OF TODD POWER SERIES 



In this section we give a proof of Theorem ^. To do this we will use Proposition ^ 
which can be used to express the coefficients fij as a sum of rational numbers, which turn 
out to equal products of certain values of the Benoulli polynomials. (Note, in contrast, that 
Proposition expresses this same number in terms of roots of unity instead of rational 
numbers.) 

Let p and q be as in the statement of Theorem ^, and let a be the cone ((1, 0), {p, q)) in 
Z^. We let pi = (1, 0), and p2 = {p, q) denote the generators of this cone. The left hand side 
of Theorem ^, fij{p,q), equals the coefficient of x^y^ in the power series ta{x,y). We now 
compute this power series using Proposition ^.3[ This proposition contains an expansion for 



the power series s^, which is equivalent to the following expansion of t^: 



U{x,y) = qxy ^ e~ 



One finds also that 



-{{■m,pi)x+{m,p2)y) 
med-nM 

Every point of a fl M can be written uniquely as a nonnegative integral combination 
of the generators ui = (g, —p) and U2 = (0, 1) of &, plus a lattice point in the semiopen 
parallelepiped 

P = {cui + du2\c, d e [0, 1)}. 
Using (mi,pi) = {u2, P2) = q, it follows that 

Ux,y) = qxy ^ V e-(<'"'''i>"+<'"'''2>^). 

^ ^ " I — Q-qx I _ Q-qy Z-^ 

mePnM 

P=l{k,{^}~^):k = 0,...,q-l\, 

[ q q J 

where {x} G [0, 1) denotes the fractional part of x. (Note this is slightly different from 
(x) G (0,1], which appears in the definition of Sij: by definition, (0) = 1, while {0} = 0.) 
One then finds that 

1 1 k kp 

k=0 

We may then compute fij{p,q) as the coefficient of x'^y^ in the above expression. It is 
convenient to replace x and y with —x and —y, which introduces a factor of (—I)*"'"-'. We 
obtain 

f],{p,q) = (-l)^+Y+^-i^coeff(x^^^ coefffy^^^ . 
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It is then clear that we can write the above sum in terms of values of the Bernoulli 
polynomials, as follows: 

= (-l)'+^g'+^-^xiBer/^)Ber,({^}). 

fe=o \ y / 

Using the identity 

Ber,(A) = (-iyBer,(l-A), 
we may rewrite our expression as 

= (-l)Y+^-^xiBer/^)Ber/(-^)]. 

The sum on the right hand side is easily seen to equal the sum defining Sij{p, q), except 
for a possible discrepancy in the k — term. It follows easily from the definitions that the 
k = terms actually match unless i = j = 1, or i = 1 and j is even, or j = 1 and i is even. 
In these cases, we need the correction terms which appear in the statement of the Theorem. 
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